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04 Abstract A monotone self -mapping of the nonnegative orthant induces a monotone 

^^ discrete-time dynamical system which evolves on the same orthant. If with respect to 

d this system the origin is attractive then there must exists points whose image under 

^. the monotone map is strictly smaller than the original point, in the component-wise 

partial ordering. Here it is shown how such points can be found numerically, leading 

to a recipe to compute order intervals that are contained in the region of attraction and 

where the monotone map acts essentially as a contraction. An important application is 

the numerical verification of so-called generalized small-gain conditions that appear 

in the stability theory of large-scale systems. 
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1 Introduction 



By IR+ we denote the nonnegative real numbers, [0,oo). A class J^ function is a 
continuous function 7: M+ -^ M+ that satisfies 7(0) = and is strictly increasing. The 
p^ function 7 is of class J(fc<, C Ji^ if in addition 7 is unbounded. Note that with respect 

f^ to composition the class J(f,„ is a group and the class J(f a semi-group. Moreover, 

iy-^ sums and positive multiples of J(f functions are again J^ functions. 

(^ The nonnegative orthant K'5_ induces a partial order on M", which coincides with 

^D the component wise ordering, and we write foix,y E R", x <y if y — x E M'j_, x <yif 

^~~^ [x <y andx y^y], andx -^y if y — x E intM1_, the interior of M']_ . 
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A map T : M.'\_ -^ M'J_ is monotone ifx<y implies Tx < Ty. For any J^oo function 
p, the map D = diag(p) : M'| -^ R'\_ defined by {Dx)i — p(x,) is an example of a 
monotone map. 

We consider the following problem: 

Problem 1 Let a monotone, continuous map T : Wl_ -^ M.1_ with r(0) =0 and a real 
number r > be given. Find s* E M'j_ satisfying 

1. Ts* < s* (which impHes s* > 0), 

2. ||5*||i =r, and 

3. Ts-^s for all s £ [0, s*],s^O. D 

The existence of such an s* (for sufficiently small r > 0) is a necessary and sufficient 
condition for asymptotic stability of the origin with respect to the discrete time system 

s+ = Ts, seR'[. (1) 

It also arises as a so-called generalized small-gain condition [16]. If it is satisfied then 
the set [0,s*] is contained in the region of attraction. Moreover, in this case also the 
system with input, 

s+ = Ts + w, s^weRl, (2) 

is locally input-to-state stable fTM], and knowledge of s* yields estimates for the sets 
of admissible inputs and initial conditions which result in bounded outputs, cf. [16). 

The numerical solution of this problem is interesting in two aspects: First of all, 
it provides a numerical way to make a qualitative assertion. Secondly, this assertion 
is not only of qualitative nature (i.e., a system is stable in some sense), but also quan- 
titative in that an estimate for the region of attraction is obtained as well. 

The numerical solution of this problem is interesting for several applications: 
In the context of large-scale interconnections of nonlinear systems the solution to 
Problemfllcan assert (local) input-to-state stability of interconnections of many sys- 
tems in arbitrary interconnection topology. This in turn can be useful for formation 
control 1 17 1 or the effective implementation of decentralized model predictive con- 
trol lfT3l . Furthermore, in the same context the knowledge of s* can be used to find 
a locally Lipschitz continuous Lyapunov function for the composite large-scale sys- 
tem. Another application is in queuing theory. Here the solution to Problemfllcan be 
used to ascertain that a given switching policy stabilizes a flow switching network 
via the use of a monotone monodromy operator, cf. |6|, at least for specified range of 
initial buffer levels and bounded inflow. 

It is known that a solution to Problem [T] must exist for any r > if the origin is 
globally attractive with respect to ([T]l. In this case necessarily it holds that Ts ^ s for 
all s e M'l ,5^0, and by virtue of a topological fixed point result, for every r > 0, 
there exists an i* with 1 1 i* 1 1 1 = r satisfying Ts* ^ s* . However, even if T : K']_ — > E'[ 
is monotone, continuous, and satisfies r(0) = as well as Ts ^ s for all s € M']_,i 7^ 0, 
the origin is not necessarily globally attractive with respect to ([T]i (but it is so locally). 
For T of particular form a sufficient condition for global asymptotic stability of the 
origin with respect to ([T]) is the existence of a diagonal mapD ~ diag(id + p): K'^ — > 
R']_ with p e J^ such that for all s 7^ 0, {Do T)s ^ s, or, equivalently {T oD)s ^ s. 
In this case the system ([2]i is input-to-state stable lfT4l[T5l . 



Also known is that if r satisfies T(s(Bv) — Ts(BTv for alls, v ^ M'^, where de- 
notes component-wise maximization, then T must be of the form {Ts)i = maxj Yij{sj) 
for all i,j, where Yjj : IR+ — > M+ are nondecreasing functions. In this case T has been 
termed max-preserving. Here, the condition Ts ^ s for all i 7^ is equivalent to the 
cycle condition, which assumes that ji^t^ o yj^i^ o . . . o 'j^-^ ,;^ o ■)^-^;^ < id for all finite or- 
dered sequences (/i,/2,. . . ,4) C {1,. . . ,«}*". If this condition holds then withe denot- 
ing the vector (1, . . . , 1)-^ and? > Gone has for^(f) ■.— ma\{te,T{te), . . . ,T"^^{te)} £ 
W\_ a continuous path q : IR+ — !• IR'j_ which is unbounded and nondecreasing in every 
component and satisfies T{q{t)) < q{t). Re-parametrization yields a path q{r) satis- 
fying ||q'(r)||i = r, which can be interpreted as a parametrized "almost" solution to 
Problem[T] cf. \§\. 

In the linear case the action of T can be represented by multiplication with a 
nonnegative matrix (i.e., every component is nonnegative), which we also denote by 
T e M."^" . It is known that the following are equivalent: 

1. ri^iforall5eM';,,s>0; 

2. there exists a D = diag(id + p) with p G J^ such that {T o D)s ^ s for all s G 
M'; , .? > (equivalently, {DoT)s^s for all .? e M'| , i > 0); 

3. the spectral radius of T is less than one; 

4. there exists a unit vector s* ^0 (with respect to the 1-norm) such that Ts* <C s*, 
hence also the ray given by rs*, r £ R+ satisfies T{rs*) ^ rs* for r > and 

Ik^lli^'-; 

5. the inverse oil —T exists and is given by the nonnegative matrix Y.k>o T''- 

The existence of the vector s* is of course related to the classical Perron-Frobenius 
theory. If T is primitive then s* is just the positive Perron-Frobenius root correspond- 
ing to the maximal eigenvalue, which coincides with the spectral radius. Further ex- 
tensions of the classical Perron-Frobenius theory exist for special classes of nonlinear 
maps, in particular for homogeneous maps and for concave maps, cf. im iTTl . 

In this paper we propose the use of a homotopy method to find a point s* satisfy- 
ing Ts* ^ s*, \\s* II 1 for any given r > 0. Our method of choice is the Kl algorithm 
proposed by Eaves as a computational version of a topological fixed point the- 
orem. There are, however, more elaborate choices of related algorithms available as 
well, cf. 12] 13 for an overview. While these algorithm have the potential of admitting 
faster convergence, they tend to be more complicated to implement. One particular 
advantage of homotopy methods is that they offer global convergence: If there exists 
a point s* with H^s* || 1 ~ r with the desired properties then it will be found. This has to 
bee seen in contrast to methods based on Newton steps, which only guarantee conver- 
gence if the algorithm is started sufficiently close to s* . Moreover, Newton methods 
usually assume some level of smoothness, whereas homotopy methods only require 
continuity. 

Once the point s* has been computed, the remaining verification of property 3 
in Problem 111 is an easy task: It only needs to be checked that the sequence T'^{s*) 
converges to zero. It is then a consequence of monotonicity that property 3 must hold. 

The paper is organized as follows. In Section [2] a few facts about monotone self- 
mappings of the nonnegative orthant and their induced discrete-time systems are re- 
called. In particular, the topological fixed point theorem by Knaster, Kuratowski, and 



Mazurkiewicz is discussed. A brief and informal description of some of the underly- 
ing principles of Eaves' and other homotopy algorithms is given in Section 3] A MAT- 
LAB version of Eaves' Kl algorithm is provided in Section!?] Section B explains a 
short procedure to solve Problemfllbased on the use of a homotopy algorithm and the 
computation of one trajectory of system ([TJ. This is followed by several numerical 
examples in Section [6] 



2 Monotone maps and monotone discrete-time systems 

This section collects a few theoretical results from the literature. 

Throughout this section let T: W\_ -^ K\_ be monotone and continuous with 
r(0) = 0. Consider also the induced discrete-time systems ([TJ and Q. Denote their 
respective solutions for initial condition s^ e K\_ and, in case of (p}, input sequence 
w = {w{k) £ M'l : fc > 0}, at time fc > by <^m{k,s^) and ^/2)(^,i , w), respectively. 
If the reference to a particular system is clear from the context we omit the reference 
to the system. Observe that both systems satisfy the ordering of solutions principle: 
If s^ < v" and for all A; > 0, w{k) < u{k) (which we abbreviate by w < u), then also 
^^{k-,s^) < 0{T](/t,v'') and0j2)(/t,5°,w) < 0j2)(/t,v'',M) for alU > 0. 

We denote the sphere in M.'\_ of radius r > with respect to the 1-norm by S,- = 
{s G E1_ : \\s\\ 1 = Y,iSi = r}. Observe that Sr is an (n — 1) simplex. 

Theorem 1 Let T : W^ — > W^ be monotone and continuous with T{0) — 0. Assume 
that the origin is attractive with respect to ([T]l and denote the domain of attraction by 
3S. Then the following assertions holds: 

1. For every s G S8, s ^Q, necessarily Ts ^ s. 

2. If R (£ R, R > is such that Sr (Z ^ then for all r G {0,R] there exists a point 
s e S,-, s >> 0, satisfying Ts ^ s. 

3. The origin is stable in the sense of Lyapunov, i.e., for every e > there exists a 
5 > such that ||s||i < 5 implies \\Ts\\\ < e. 

A proof can be found in lfT4l . The first assertion is not difficult to prove directly, and 
the third follows from the second. The second assertion is the most technical, and it 
is based on the covering theorem by Knaster, Kuratowski, and Mazurkiewicz (KKM) 
llOI . Our reasoning is based on the extension given in [12 J which allows to consider 
coverings of open instead of closed sets. The argument is basically the following: 
One has Ts ^ s for all s G S,-, and Sr is a simplex. This ordering condition implies 
that the simplex is covered by the sets 

Qi = {seSr: {Ts)i<Si}. 

Moreover, no Q, can contain the face (7, C 5, opposite of the vertex re, £ f2,, where 
e, denotes the /th unit vector On the other hand, every ^-dimensional simplex a — 
co{re,},£|,| ;^i, fe < n— 1, is contained in the union [JieU, u}^'- These are the 
prerequisites of the KKM Theorem, which then ascertains that the intersection f] 12, 
must be nonempty. The proof of the KKM Theorem is based on the fact that every 
simplicial refinement of Sr must contain at least one special simplex, which is again 



covered by all sets 12, but not by any strict subclass of {i2i, . . . ,i2„}. As the size of 
the refinements tends to zero, this special simplex contracts to a point — the point of 
interest s* E S,- satisfying Ts* -^ s*.ln essence this proof technique relies on a fine 
discretization of S,- and an exhaustive search, which is not implementation friendly if 
n becomes large. 



3 Homotopy based fixed point algorithms 

An alternative proof of the KKM result has been given in fS] by means of a homotopy 
algorithm. In contrast to the original proof, here the initial simplex is successively 
refined in every iteration, thus the area that must contain s* becomes smaller and 
smaller as iterations progress. We give a simplified account on the ideas behind this 
algorithm. 

As in the original KKM paper, each point in 5 e 5,- is assigned an integer label 

l{s)^miix{i: seQi}. (3) 

Observe that /(re,) = / for / = 1, . . . ,«. To keep things uncluttered, let us call a set 
of k distinct points (vertices) a = {v' £ R"^}^^^ a {k—1 )-set. The barycentric centre 
of such a (fc — l)-set is the point v = j ^f^j v^. We call an (n — l)-set a complete if 
all vertices have distinct labels. Equivalently, every label !,...,« gets assigned to a 

in 
H=l 



vertex exactly once. Observe that the (n — l)-set ff" = {re,}JLj is complete if Ts ^ s 
for all s E Sr- 

Now one can also consider an n-set T obtained from a given (n — l)-set a by 
adjoining one additional vertex v taken from the convex hull co(c7) of the vertices of 
<T. For example, one could augment a with its barycentric centre to obtain such an 
n-set T. 

The following observation is at the core of the fixed point algorithms by Eaves 
and also at the core of aU related homotopy algorithms, cf also [2 Thm. 1.7]. 

Tlieorem 2 Every n-set has either none or exactly two complete [n — \)-subsets. 

Proof Let T = { v^ , . . . , y"+' } denote the given n-set with its n + 1 vertices. Now either 
/ ({y' , . . . , y"+' }) contains the set { 1 . . . . , n}. In this case one label must get assigned 
twice, i.e., there exists a unique / and j^k such that l{v^) = l{v'') = i. In this case 
T\{y-'} and T\y'^ are both complete (n — l)-sets. All other (n — l)-subsets must 
contain y^ and y^ and hence cannot be complete. And if {1 , . . . , n} ^ Z ({y' , . . . , y"+' }) 
then no such (n — 1 )-set can exist. D 

The basic algorithm is now the following: ff" = {'"^i}"=i serves as a complete 
entry (n — l)-set. Successively, another vertex is added, say the barycentric centre y 
of a", to obtain the n-set t — a^lJ {v}. By Theoremol T contains exactly one (n — 1 )- 
subset distinct from a^, which we denote by a' . Progressing inductively, one obtains 
a sequence of complete (n — l)-sets ff^' whose area tends to zero. 

The catch, however, is that the sequence does not necessarily contract to a point. 
Instead, a'^ may become "long and thin" as k becomes large. To prevent this kind 



of behaviour, a more sophisticated choice of new vertices v is required. For this 
choice there exist a variety of alternatives. The paper by Eaves |5| proposes two 
such choices, which are named Kl and K2. Both guarantee convergence. Of these 
two choices Kl is the easiest and shortest to implement, and therefore we have cho- 
sen Kl for this exposition. On the other hand, it does not converge as quickly as K2, 
as has akeady been observed in |5|. It should be noted, however, that even more so- 
phisticated pivoting strategies and restart algorithms can be found in |I2B1 and the 
references contained therein. A detailed discussion of these is far beyond the scope 
of this paper. In principle any of these could have been used instead of our particular 
choice for Eaves' Kl algorithm here. 



4 Implementation of Eaves' Kl algorithm 

The integer labeling function (|3]l is numerically not feasible, instead we use the func- 
tion 

le{s)=miix{i:iTs)i + e<Si}, (4) 

where e > is a design parameter. It guarantees that if a point s* is obtained with 
the fixed point algorithm, then this point does in fact satisfy Ts* <C s* and not just 
Ts* <s*. 

Obviously, if a point s* with Ts* <C s* exists at all, then it has to be found by the 
algorithm if only e > is small enough. On the other hand, from practice it is fair to 
say that larger e yield faster convergence (i.e., fewer iterations are necessary to obtain 

n- 

Also it could be noted that instead of the maximal / in Q the minimal or even any 
other unique choice should in theory do equally well. In particular, this may give rise 
to a different pivoting strategy by shifting the pivoting from the homotopy algorithm 
to the labeling function. 



4.1 MATLABcode 



Listing 1 Eaves' algorithm based on the Kl complex fSl implemented in MATLAB. 

1 function [kkmpt , noit , succ] = ea ves Kl (monmap, r, n) 

2 succ = 0; noit = Inf ; global eavesKI MAXREFINE 

3 if \sempiy (eavesKI .MAXREFINE) , eavesKI _MAXREFINE=^ e3 \ end 

4% Initialization of tine complex 

5 q = -eye{n) + diag( ones (n-1 ,1) , -1); i/ = zeros (n,n + 1 ); 

6 1/(1 ,1) = 1 ; v[n,^) = ^ ; g=[n,1 :n-1]; 

7 for / = 1 : n, i'{:,/ + 1) = v(: , i) + q{:,g{i)); end 
s% compute initial labels 

9 istar = ^■, inew=^^, /=zeros(1 ,n+1 ); 

10 for / = 1 : n+1 , 

11 X = v{: , i )/ norm ( i/ (:,/), 1 ) * r; 

12 try 

13 '(') = eavesLabel (monmap, x) ; 

14 catchi exception , 

15 error (' initial^label^computation^went^wrong' ); 



16 end 

17 end 

18% the main loop 

19 lexo = ^ ^.n+^ ; 

20 for noit = ^■.eavesK1.MAXREFINE, 

21% identify vertex to drop and add new one 

22 f=find{ /==/(/new)); 

23 if f(1) == inew, istar=t{2); else istar^t (^)■, end 

24 if istar==^, vnew=v{: , lexo{n+^)) + q{^. ,g{^))\ inew=n+^^, 

25 lexo = [lexo{2^.n+^) lexo{^)]^, v{: , lexo{inew)) = vnew; 
g=[g{2:n) g(1)]; / = [/(2:n+1) 0]; 

27 elself istar < n+^ , 

28 vnew=v{: , lexo{ istar -^)) + q(■. ,g( istar)); 

29 inew=istar; lexo=lexo; v{: , lexo{inew)) = vnew\ 

30 temp = g{ istar); g( istar) = g( istar -^); 

31 g{istar -^)=temp; clear temp; 

32 / = [/(1 : /sfar-1) I { istar + ^■.n + ^)]; 

33 else i/neiv=i/( : , /exo{1 )) -q{ : , g(n) ) ; /neiv=1; 

34 lexo = [lexo{n+^) lexo{^■.n)]; v{: , lexo{inew)) = vnew; 
g=[g{n) g(1:n-1)]; / = [0 /(1:n)]; 

36 end 

37% compute label of new vertex: 

38 l(l<mpt = vnew/ norm( vneiv, 1 ) * r; 

39 try 

40 I {inew) = ea\/esLabe\ {monmap, kl<mpt); 

41 catcli exception , 

42 error { ' successive^label^computation^went^wrong' ) 

43 end 

44 if all { feval (monmap, /</cmpf)</(/<mpf) , 

45 succ = 1 ; % success I 

46 return ; 

47 end 

48 end % refinements 

49 end % function eavesKI 

50 function / = eavesLabel (monmap, x) 

51 global eavesLabel-distance 

52 if \sempty(eavesLabel. distance) , eavesLabei. distance = 1e-2; 

53 end 

54 Tx = f eval (monmap, x) ; 

55 f=find(rx+ eavesLabel-distance<^x); /=max(f); 

56 if isempty(/), error ([ 'Found^point^without^label.^Try^' ... 

57 ' decreasing^eavesLabel_distance . ' ]); end 

58 end % function eavesLabei 



Of the output arguments kl^mpt denotes the vector s* £ M'|, noit the number of 
iterations that have been consumed by the algorithm to find l<kmpt, and succ is either 
1 to denote that the algorithm was successful and kl^mpt is a point of interest and 
otherwise. 

Of the input arguments monmap is a function handle to the monotone map 
T: W\_ -^ M!\_ which satisfies T{Q) = 0. The parameter r is the radius r > of the 
sphere Sr C M!]_, where the algorithm tries to find s* . Lastly, n denotes the dimension 

ofM';. 

There are two additional global variables that can be tweaked to modify the per- 
formance of the implementation. They are eavesKI .MAXREFINEwith a default value of 



1000 and eavesLabeLdistance with a default value of 10 ^. Their meaning is explained 
below. 

4.2 Usage 

To use the above MATLAB function, one has to implement the monotone map T of 
interest into a MATLAB function. An example of this is given in Listing [2] 

Listing 2 A MATLAB implementation of the monotone map defined in Examplelll cf. Sectionl6| 

1 function y = T{x) 

2 y=zeros(size{x)) ; 

3 for /<=1 :length (x)-1 , 

4 y(k) = y{k) + x(k+^)'{k+^)■, 

5 y{k+^) = y(k+:) + x(/<)"(i /(/t+i )); 

6 end 

7 y=y/4; 

8 end 

Now, to compute a point s* e R1_ satisfying Ts* ^ i*, ||i* || i = r > 0, one calls 

[s.star , noit ,succ] = eavesKI {@7; /-,n) 

which should yield the desired result. If necessary, the behaviour of the implementa- 
tion can be fine-tuned by specifying 

global eavesKI .MAXREFINE eavesLabeLdistance 

and assigning new values to these variables for the maximal number of iterations, and 
respectively, the parameter e appearing in the labeling function Q. 



4.3 Convergence 

Assuming arbitrary precision computations and number representation as well as suit- 
able choices for eavesKI .MAXREFINE (maximal number of iterations) and 
eavesLabeLdistance (i.e., e > in (|4]i), the algorithm does what is expected: 

Theorem 3 Let T : M'J_ — > M.']^_ be monotone and continuous. Let r > and assume 
that Ts ^ sfor all s G 5,-. Assume that there exists a point s* € Sr such that Ts* <C s*, 
< e < minjij" — {Ts*)i}, where e is the parameter in Q. Then the algorithm in 
ListingU\produces a point s** G Sr (possibly different from s* ) satisfying Ts** ^ s**, 
provided that the maximal number of allowed iterations is large enough. 

Proof The claim follows from the corresponding more general result in ||5l. 



4.4 Remarks on computational complexity 

In each iteration of the algorithm, the map T : W\_ -^ W\_ has to be evaluated once to 
compute label of the new vertex. In addition, n comparisons are necessary to find the 
old vertex with the same label as the new one to be dropped. 



5 Algorithmic solution to Problem [T] 

Building upon an implementation as in the previous section. Problem [T] can now be 
solved as follows: 

1 . Given r > compute s* using the algorithm in the previous section, or using a 
more sophisticated implementation based on one of the algorithms proposed in 
e.g. HJ-HJ. If this is successful then properties 1 and 2 of Problem [T] are already 
satisfied. 

2. Compute {T^s*}i;>(). If this is a null-sequence then Ts ^ s for all s e [0,5*], 5 ^ 0, 
i.e., property 3 of Problem [T] is satisfied. 

The assertion of the second step in this short meta-algorithm relies on the first 
statement in Theorem n] For if T''s* tends to zero as ^ ^- then s* e ^, the region 
of attraction. By the ordering of solutions principle also every point s e [0,s*] must 
belong to 3§ as well. 

The computational complexity of the second step consists of the computation 
of only one trajectory of a discrete-time monotone system. The trajectory has to be 
bounded, because due to the ordering of solutions principle it must be confined to 
the order interval [0,^*]. Furthermore, due to the monotonicity of T it has to be non- 
increasing in every component, which allows to terminate further computation once 
0m (A;, i*) is sufficiently small for some k>Q. 



6 Examples 

In this section we consider a few numerical examples. The first is a nonlinear map 
T: M'5_ -^ M'l which can be defined for any « > 2. For this map T it is known that 
system (|2| is input-to-state-stable, implying that the origin is globally asymptotically 
stable for system ([T]), and hence the Eaves Kl algorithm should produce an s* for 
arbitrary n > 2 and r > 0. 

The second example is a statistic generated from randomly chosen nonnegative 
matrices A E W^" with spectral radius less than one. 

A third example shows that for a given x{Q) the pure iteration ofx{k+ 1) = Tx(k) 
does in general not produce an x{k) such that Tx{k) <^ x{k), even if fc > is large. 

Example 1 Let n > 2. Consider the nonlinear map T : M'| -^ M']_ defined in IfTSl Ex- 
ample IV. 1 ] given by 

with the convention that sq — i^+i — 0. For example, in the case n = 5 one has 
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Observe that T{Q) = 0, and that T is obviously monotone and continuous. It has 
been shown in fTS^ Example IV. 1] that the induced system (|2]i is input-to-state stable. 
Moreover, it can be verified directly that for every r > 0, the vector 



p{r) 



KrrJ 



satisfies Tp{r) <C p{r)- Since p: IR+ -^ M!\_ is continuous, and in every component 
unbounded, for every r > there must exists a point s* G 5, such that Ts* ^ s*. 
Figure [T] shows how long it takes (in terms of iterations) for the Eaves' algorithm to 
find such a point 5*. 
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Fig. 1 Number of iterations needed by Eaves' Kl algorithm for r = 10, and different choices of n and 
e > in ExamplehJ The maximal number of iterations was set to 10^. It took 67 seconds to generate this 
whole plot on an 2.4 GHz Intel MacBook. 



Example 2 It is relatively easy to generate many positive nxn matrices with a spec- 
ified spectral radius in MATLAB. If the spectral radius p (A) of a nonnegative matrix 
A is less than one then it defines a monotone mapping that should allow for a point 
s* S Sr with Ai* ^ s*. Here we have chosen the spectral radius to be p{A) = 0.8 and 
have generated a number of matrices A G W^" for different choices of n. Again we 
have applied the algorithm in Listing[T]and counted the number of iterations needed. 
The outcome is plotted in Figure l2] 



Example 3 Consider the monotone map T : 



T{x) 



x2,Xx]) , with X e (0,1). Obviously T{0) 



— > M^ given by 
and it is easy to check 
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Fig. 2 Mean number of iterations needed for linear maps with spectral radius p{A) = 0.8 together with 
minimal and maximal number of iterations needed. Each data point corresponds to 10 randomly chosen 
matrices (using MATLAB's rand(n)). Generating this plot took 205 seconds (not taking into account 
generation of the test matrices). The maximal number of iterations was set to lO'. 



that for any x S W\_, T'^x ^> as A: — )> oo. Hence, T is a contraction. Yet, given x G M^, 
if not already Tx-^x then there exists no ^ > 1 such that r*+'x <C T'^x. For we have. 



T^x^ 



/Xx\ 

Xx2 



Tx = 



X2 



Assuming that Tx <^ x, we have either {Tx)i — y/x2 > x\, which impHes (r^x)2X 



Xx2 > Axj = {Tx)2. Otherwise, we have {Tx)2 — Xx\ > X2, implying that (T x) 



vAxi > y/x2 — {Tx)i. This deduction repeats inductively. As a consequence, a pure 
iteration of the map T cannot yield a solution to Problem [T] 



7 Conclusions 
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